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Abstract —We address the issue of applying quantized com¬ 
pressed sensing (CS) on low-energy telemonitoring. So far, 
few works studied this problem in applications where signals 
were only approximately sparse. We propose a two-stage data 
compressor based on quantized CS, where signals are compressed 
by compressed sensing and then the compressed measurements 
are quantized with only 2 bits per measurement. This compressor 
can greatly reduce the transmission bit-budget. To recover signals 
from underdetermined, quantized measurements, we develop a 
Bayesian De-quantization algorithm. It can exploit both the 
model of quantization errors and the correlated structure of 
physiological signals to Improve the quality of recovery. The 
proposed data compressor and the recovery algorithm are vali¬ 
dated on a dataset recorded on 12 subjects during fast running. 
Experiment results showed that an averaged 2.596 beat per 
minute (BPM) estimation error was achieved by jointly using 
compressed sensing with 50% compression ratio and a 2-bit 
quantizer. The results imply that we can effectively transmit n 
bits instead of n samples, which is a substantial improvement 
for low-energy wireless telemonitoring. 

Index Terms —Quantized Compressed Sensing, Block Sparse 
Bayesian Learning, Data Compression, Telemonitoring 


I. Introduction 

In wireless health monitoring systems [1], large amount and 
various types of physiological signals are collected from on- 
body sensors, and then transmitted to nearby smart-phones or 
data central via wireless networks. Energy consumption is a 
critical issue in these systems. Compressed sensing (CS) [2] 
is a promising technique for such systems for low-energy data 
acquisition, compression, and wireless transmission [3], [4]. 

In the framework of CS, a signal x G is compressed 
by a simple matrix-vector multiplication, 

y = + n. (1) 

where # £ j^mxn called the sensing matrix, y £ 
is the compressed measurements and n is the measurement 
noise. Usually is underdetermined, i.e., M < N, and the 
ratio (iV — M)/N is called the compression ratio of CS. 
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Compression via CS can be implemented in two ways as 
shown in Fig. 1. One is analog CS [5]-[8], where the matrix- 
vector multiplication $x is implemented in analog domain and 
the compressed measurements are quantized via an Analog-to- 
Digital Converter (ADC). Analog CS is also called Analog- 
to-Information Converter (AIC) [7], [8]. Usually, an AIC is 
implemented in a dedicate ASIC chip to achieve low-energy. 
The other is digital CS, where signals are firstly quantized into 
digitalized samples and $x is calculated in MCU [9] or FPGA 
[10]. The advantage of digital CS is that it can utilize exiting 
ADCs and can be implemented as an efficient, low-energy data 
compressor [10]. In both cases, the sensing matrix $ needs 



Fig. 1. Compression via Compressed Sensing (CS). The matrix-vector 
multiplication in (1) can be either implemented using analog circuit in (a) 
or digitally in (b). In both cases, an ADC is required to convert continuous 
signals to quantized numbers. 

not to be transmitted and is stored in both the transmitter and 
receiver. 

At a receiver, the original signal x is recovered from the 
compressed measurements y via 

min||x||o subject to ||y—4>x| I2 < £ (2) 

where ||x|jo is the norm of ||x||, and e is the tolerance 
of noise or modeling errors. Calculating the solution is very 
hard. Generally, one seeks the solution of a relaxed convex 
optimization problem [11], in which ||x||o is replaced with 
||x||i or other terms encouraging sparse solutions. The quality 
of signal recovery can be improved by exploiting other struc¬ 
ture information in signals, such as wavelet-tree structure [12], 
piecewise smooth [13] or block sparse structure [3], [14]. 

CS has been successfully used in low-energy telemoni¬ 
toring of physiological signals such as EEG [15] or fetal 
ECG [3]. However, these works [3], [15] assumed that the 
compressed measurements y was real-valued. But in practice, 
the compressed measurements y must be quantized before 
transmission, i.e., 

y = Q($x) 


(3) 
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where Q(-) is the quantization operator that maps a real¬ 
valued signal to finite quantization levels [16]. This process un¬ 
avoidably introduces errors, called quantization errors. Some 
CS algorithms [17]-[19] were developed to recover signals 
from the quantized measurements. This recovery procedure is 
often called de-quantization [20] and a CS algorithm is also 
called a decoder [21]. Haboba et al. studied the quantization 
errors and its effect on signal recovery using synthetic signals 
with fixed sparsity [7]. However, the signal model used in 
[7] is not practical for telemonitoring applications, as most 
physiological signals are only approximately sparse. Wang et 
al. analyzed the performance-to-energy trade-offs introduced 
by quantized CS [8] in EEG telemonitoring and provided a 
brute force method searching the optimum configuration of the 
quantizer. But the CS algorithm used in [8] ignored the model 
of quantization errors during signal recovery, which may yield 
degraded performance. 

Although quantization introduces errors, suitably using 
quantization can largely reduce the wireless transmission bit- 
budget [22]. In this work, we address this issue and use 
quantized CS for low-energy telemonitoring. Eirst, physio¬ 
logical signals are compressed according to (1), and then 
quantized according to (3) with 2 bits per measurement. 
This compression scheme greatly reduces the transmission bit- 
budget, which benefits to low-energy telemonitoring. 

On the de-compression and de-quantization stage, we pro¬ 
pose a Bayesian de-quantization algorithm, denoted by BDQ. 
It exploits correlation structure within physiological signals 
and also takes into account the quantization errors. Note that 
this algorithm does not exploit sparsity to recover signals, as 
done by most other compressed sensing algorithms. Instead, it 
exploits correlation of physiological signals. Our motivation is 
that during wireless health monitoring many raw physiological 
signals are usually less sparse, namely these noisy signals are 
not sparse in the time domain and also not sparse in other 
transform domains [23]. In this case exploiting sparsity may 
not be very effective. In contrast, exploiting correlation may 
be a better direction, as shown in [3]. However, the work in [3] 
did not consider the quantization errors, and thus has inferior 
performance to our proposed BDQ algorithm, as shown in 
experiments later. 

We study the application of Photoplethysmography (PPG) 
and accelerometer telemonitoring for fitness training, in which 
the heart rate must be accurately estimated during intensive 
physical exercises. We exploit the optimum compression by 
jointly tunning the compression ratio (CR) of CS and the 
quantization bit-depth. The experiment results show that an 
averaged 2.596 BPM absolute heart rate estimation error is 
achieved with CR = 0.50 and a 2 bits quantizer. The Pearson 
correlation between the estimated heart rate from the recovered 
datasets and the ground-truth is 0.9899. These results imply 
that we can effectively compress raw segments of PPG and 
accelerometer data from N samples into N bits for low-energy 
telemonitoring. 

The rest of the paper is organized as follows. Section II 
introduces the quantization model of compressed sensing. The 
Bayesian De-Quantize algorithm are presented and discussed 
in Section III. Section IV describes the experimental set up 


and numerical results. Discussions are given in Section V and 
Section VI concludes the paper. 

II. The Quantization Model oe Compressed 
Sensing 

The function of an ADC is time sampling and scale quanti¬ 
zation, which is illustrated in Pig. 2. Let B denotes the number 



Fig. 2. Time sampling and scale quantization. The quantization error is 
denoted by and the cell width is A. The mid-point value within a cell is 
taken as the quantized value for a sample point falling in that cell. 


of bits per measurement, which is also called the bit-depth 
of a quantizer. Represented by B bits, the scales between 
the positive and negative reference voltage [—14-e/, 14-e/] are 
divided into L = 2^ quantization levels. The cell width [16] 
A of a uniform quantizer is 


_ 2Vref 

2B 


(4) 


Signals larger than Vref or smaller than —Vref are saturated. 
Those saturation signals are quantized with the same level as 
signals lie in [Vref - A, Vref] or [-Vref + A, -Ere/]- 

Por a sample S K falls in a cell, the mid-point value 
in that cell is used as the quantization value Q(vi). The 
quantization error is defined as 


Ci = Qiv,) - Vi, 


(5) 


Ci distributes uniformly between [—A/2, A/2]. Let Ve = 
[—A/2, A/2] denotes the domain of quantization errors, then 


p{ei)^U{Ve). 


(6) 


It should be noted that Ve is unbounded when a signal 
saturates [19]. Purthermore, the variance of the quantization 
error, denoted by (Je, is. 


tTe = 


12 ■ 


(7) 


A. Analog CS 

Por analog CS, the matrix-vector multiplication $x is 
implemented in analog domain and the compressed measure¬ 
ments are quantized via an ADC before transmission, 

y = -f n, y, = Q(y) 


( 8 ) 
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where the subscript q denote the quantized signals. Let e = 
Yq - y, we have, 

Yq = f>x + n + e. (9) 

The transmission bit-budget of analog CS is MB bits, which 
can be controlled by the number of compressed measurements 
M as well as the quantization bit-depth B of the ADC. 


B. Digital CS 

For digital CS, analog signals are firstly quantized by an 
ADC via Nyquist sampling, then the digitalized signals are 
compressed in MCU [9] or FPGA [10] via a matrix-vector 
multiplication, 

y = $x,-fn, x^ = Q(x) ( 10 ) 

In practice, larger bit-depth Bi is required for the ADC to 
reduce the distortions during signal acquisition (i.e., Bi = 12 
was used in [10]). For larger Bi, the variance of the quanti¬ 
zation error e is much smaller and ( 10 ) is simplified as, 

Y^^Xq + n, X, «x (11) 


C. The two-stage data compressor 

In (10), both y and Xq are represented in fixed-point arith¬ 
metic with Bi bits [10]. We may further reduce transmission 
bit-budget by simply rounding each sample of y to i? bits, 
denoted by Yg = Qr(y), where Qr(-) is called a rounding 


y 

a = i 
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Fig. 3. The operator Qr(-) can be efficiently implemented in digital CS. In 
this example, the compressed measurements y are represented using unsigned 
fix-point arithmetic (Bi = 8 and the fraction length is also 8 bits). The 
rounding outputs have only 2 bits per sample. 


operator and its function is illustrated in Fig. 3. Qr can be 
regarded as an economy quantizer in fixed-point arithmetic. 

We denote by = Yg ~ V the rounding error, which 
assumed to be uniform distributed. 


( 12 ) 

where = [—Ar/ 2 , Ar/ 2 ] and the cell width A^ = 

After rounding, (11) can be reformulated as. 


Yg=^Xq-(-n-\-er, X, « X (13) 


The bit-compression ratio CRt, i.e., the ratio between the 
reduced bit-budget after compression divided by the total input 
bit-budget, is defined as 


CRf, 


NBi - MB 
NB, 


1 - (1 - CR) • 


(14) 


where CR = {N — M)/N is the compression ratio in CS- 
based telemonitoring [3], [10]. 



Fig. 4. A unified framework for low-energy wireless telemonitoring using 
quantized compressed sensing. 


In our paper thereafter, we do not distinguish between 
quantization Q(-) in analog CS (9) and rounding Qr{-) in 
digital CS (13). Instead, we model them in a unified framework 
in Fig. 4, which is a two-stage data compressor that can be 
formulated as, 

z = $x -f e -f n (15) 

where z is the quantized measurements represented with B 
bits per sample, $ S sensing matrix, n and 

e are measurement noise and quantization noise respectively. 
In this framework, we transmit only MB bits instead of N 
samples for wireless telemonitoring. 

III. The Bayesian De-Quantize Algorithm 
A. Bayesian Hierarchical Model 

1) Noise Model: The measurement noise n is usually 
assumed Gaussian with variance A, i.e., n ~ AI). From 
(15) we have 

p(z|x, e; A) = A/'(4*x 4 - e, AI) (16) 

The quantization error e is uniform distributed, 

p{e)=U{Vg) (17) 

where Vg = [—A/2, A/2] and A is the cell width of a 
quantizer. In practice, the cell width A is known a prior given 
B and the reference voltage Lre/- 

Remark 1: The quantization error models rounding errors 
between the analog input and the digitalized output. It is non¬ 
linear especially for low-resolution ADCs and also depends 
on the amplitude and frequency of a signal. Exploiting the 
dependencies between signals and quantization errors may 
improve the quality of recovery, however it is difficult to 
do so [19]. To simplify our model, we assume an uniform 
distributions for e and do not consider the dependency between 
the quantization error and the analog input. 

Remark 2: We studied only multi-bit (B > 2) quantized 
CS and do not consider the extreme case of 1-bit compressed 
sensing. 1-bit CS may have better recovery performances than 
multi-bit CS under the same compressed bit-budget MB [18]. 
Flowever, 1-bit CS loses scale information of signals, which 
challenges its use in low-energy telemonitoring applications. 
We refer the reader to the literature on 1-bit compressed 
sensing [18] for more details. 

2) Signal Model: We assume a correlated structure within 
the signal [14] where we model the prior of signal x as, 

p(x| 7 ,P) =AA( 0 , 7 P) (18) 

where 7 is a non-negative parameter controlling the variance 
of the signal x, P is a symmetric positive semi-definite matrix 
modeling the correlation structure of signal x. The diagonal 
entries of P are normalized to Is during iterative learning. 
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B. The Bayesian De-Quantization Algorithm 

We estimate {x, e, A, 7 , P} using their joint MAP estimator, 

{x,e,A, 7 ,P} = argmaxlogp(x, e, A, 7 ,P|z) 

= argmaxlogp(e|z,x, A, 7 , P) • p(x, A, 7 , P|z) 

A nested Expectation Maximization (nest-EM) approach 
[24], [25] is adopted. The nest-EM is an iterative monotoni- 
cally convergent method, it consists of an inner and outer EM 
loop, which is briefly sketched in Eig. 5. 



Fig. 5. Bayesian De-Quantization Algorithm using Nested EM approach. 


We initialized the nest-EM procedure with eg = 0 and let 
y = z — eg, then iterative over, 

• Inner E-Step, estimate and Sj, from posterior prob¬ 
ability p(x|y; 7 ,P, A), 

• Inner M-Step, update 7 , P and A via maximize the 
likelihood ^( 7 , P, A|y), 

• Gutter E-Step, estimate the first moment of quantization 
error e from p(e|z, x; A), and update y = z — e. 

1) Inner-E step: The posterior p(x|y; A, 7 , P) can be ex¬ 
pressed in an analytical form [26], 


p(x|y; A, 7 ,P) = 


p(y|x; A)p(x| 7 ,P) 

p(y|A,7,P) 


= S-'(x - rtj , 


where p, and S are respectively. 


/x = 7P$^(AI + 7^P$^)-V (19) 

S = (-P-i -f (20) 

7 ^ 


2) Inner-M step: The parameters A, 7 and P are estimated 
by a Type II maximum likelihood procedure [26], 


/:(A,7,P) = -logp(y|A,7,P) (21) 

= log |AI + 7 $P$^| -f y"^ (AI + 7 $P$^)”^ y (22) 


Optimize over £(A, 7 , P), we have update rules for A, 7 and 
P respectively. 


^ ||y-$/x||i + Tr(S$^$) 

M 

(23) 

p _ S + pp^ 

(24) 

1 


7=^Tr[p-i(S + /r/x7] 

(25) 


Remark 3: Due to the coupling of e and n in (15), the 
estimate of the noise variance A is inaccurate due to an 
identifiability issue as stated in [27]. Therefore, A is often 
treated as a regularize parameter [10], [19]. We set A = 0.001 
as the default parameter used in the experiment. 


3) Regularization on P.- Regularization on P is required 
due to limited data [14]. In [14], the author provided an 
empirical method on the regularization of P using a symmetric 
Toeplitz matrix, 

P„=rl*-^l, = |r|<l (26) 

where r is the correlation coefficient empirically calculated 
from the ratio between the mean of sub-diagonal of P and the 
mean of main diagonal of P. Such regularization is equivalent 
to modeling the correlation structure as a first-order Auto- 
Regressive (AR) process [27]. 

An AR(1) matrix has simple tri-diagonal inverse, 

/I -r 0 ••• 0 \ 

—r 1 -I- —r • • • 0 

P"' = r^ ■ ■■■ ■ ^27) 

1 — ' 

0 ••• —r l-fr^ —r 

\ 0 ••• 0 -r 1 / 

and it can be decomposed as 

P“i = -^T (D^D) (28) 

where T = diag{l, • • • , 1,1/(1 -f r^)} is a diagonal matrix 
and D is a temporal smooth operator [28] defined as 

/I -r 0 ••• 0\ 

0 1 -r • •• 0 

D- , (29) 

0 0 ••• 1 -r 

yo 0 ••• 0 1/ 

A similar regularization method was proposed in [29] which 
models the inverse covariance matrix as D^D. 

This type of regularization can exploit the correlated struc¬ 
tures within physiological signals [3], [15]. However, the 
correlation coefficient r in (26) or (29) can only be empirically 
calculated [14] or fixed [29] in existing algorithms. 

We presented in this paper an AR(1) approximation to 
estimate the correlation matrix in (24) using Karhunen-Loeve 
Transform (KLT) [30]. The correlation matrix P calculated by 
(24) is a symmetric, positive semi-definite matrix. Let 

UDU"^ = SVD(P) (30) 

denotes the Singular Value Decomposition (SVD) of P, where 
columns of U are the eigenvectors of P and D is a diagonal 
matrix with d = diag(D) being the descending ordered 
eigenvalues of P. It have been shown in [31] that the basis 
vectors of a Discrete Cosine Transform (DCT) approach to 
the eigenvectors of the inverse of an AR(1) matrix as the 
coefficient r goes to 1. More precisely, from (27) we define 

K^P“i(l-r2), (31) 

then the rows of a DCT Type-2 matrix is the eigenvectors 
of K as r —)• 1. This property of DCT has made it a pop¬ 
ular transform for decomposition of highly correlated signal 
sources [30]. 
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From (24), (30) and (31), we may regularize the correlation 
matrix P with an AR(1) matrix by simply substituting the 
eigenvalue matrix U in (30) with a DCT Type-2 matrix, 

P = dctmtx(iV)^ D dctmtx(iV) (32) 


where dctmtx(A^) generates a DCT Type-2 matrix of size N. 
We then normalize the diagonal entries of P to Is by 


^ = \J diag(P), P = diag ^(v)Pdiag ^(v) (33) 

where diag“^(v) builds a diagonal matrix with the diagonal 
entries given by v. In (33), P was the AR(1) approximation 
to the correlation matrix P. 

7 was updated afterwards using the regularized P, 

7= ^Tr[p-i(S-l-/x/r^)] (34) 

Remark 4: Using KLT to exploit the correlation structure 
in signals has long been exist in literature. In [30], the author 
studied the use of Toeplitz matrix to approximate empirical 
correlation matrix. To our best knowledge, our work was the 
first to use KLT in regularized least squares to solve under¬ 
determined optimization problems. The DCT approximation 
in (32) is also attractive for its computational efficiency in 
calculating /r, S and 7 in (19), (20) and (34). 

4} Estimate The Quantization Error e (Outer-E Step): We 
calculate the expected value of e from 

.A/2 

e = ^[ele G I?e] = / e p(e|z, x; A)de 
J-A/2 


which is the expected value of a truncated Normal distribution, 
e can be obtained analytically [19], 


G — /Xg 


^ PDF(le) - PDF(ug) 
CDF(lg) - CDF(ug) 


(35) 


where /Xg = z — le = (~A/2 — fie)/V~\ and Ug = 

(A/2 — /Xg)/-\/A. PDF(-) and CDF(-) are the probability 
density function (PDF) and cumulative density function (CDF) 
of a standard Normal distribution respectively. 

Remark 5: Besides the quantization cell width A, we may 
also have the prior information of the reference voltage lAe/- 
For the estimates of the measurements z = values larger 
than Vref might indicate that saturation occurs. In this case, 
we may write (35) as. 


Gl — /Xgl — a/A • 


PDF(lgT) 
CDF(lgx) - 1 


(36) 


where I denotes the set of index where z > Vref- 

5) The proposed algorithm: The resulting algorithm is 
summarized in Fig. 6, named as the Bayesian De-Quantize 
algorithm (BDQ). 

Remark 6: The BDQ algorithm can be used to recover 
piecewise smooth signals from quantized, and possibly under¬ 
determined measurements. It shares some similarities with the 
Block Sparse Bayesian Learning (BSBL) framework [14] in 
non-sparse mode where only the correlations within signals 
are exploited. However, in the experiment we found that the 
regularization method (32)-(33) on P in BDQ was superior to 
the empirical methods in BSBL, which yielded better recovery 
results on physiological signals. 


1 : procedure BDQ(z,$) 

2: Outputs; X, e 

3: Initialize: 7 = 1, P = I, A = 0.001 

4: while not converged do 

5: Estimate /x, S by (19), (20). 

6: Calculate P via (24) 

7: Regularize P via (32)-(33). 

8: Calculate 7 via (34). 

9: Calculate e by (35) 

10: Update y = z — e 

11 : end while 

12 : end procedure 

Fig. 6. The Bayesian De-Quantize (BDQ) algorithm. 


IV. Experiments and Results 

A. Datasets 

We simultaneously collected ECG, PPG and accelerometer 
signals from 12 volunteers with age ranged from 18 to 35. 
Pig. 7 shows the hardware set up for data recording. Por 



Fig. 7. The hardware setup for data recording. ECG, PPG and accelerometer 
signals were simultaneously collected. Wet electrodes were used to obtain a 
three-lead ECG, where the leads LA and RA were placed at left and right 
chest respectively, and the LL lead was placed at the left lower abdomen. 
PPG signals were collected using reflective pulse oximeter sensors with green 
LED (wavelength: 515nm). Two PPG sensors were placed at the back of the 
wrist while another one was placed right on the pulse position. Three axial 
accelerometer data were collected using ADXL345 with ±4g range. Data 
were transmitted to a laptop via bluetooth. 


each subject, we recorded data for 5 ~ 6 minutes. During 
data collection, a subject ran on a treadmill with speeds 
ranged from 2km/hour to 15km/hour. The sensor band was 
conveniently worn on the wrist and we intentionally introduced 
additional artifacts by asking all subjects to pull clothes, wipe 
sweat, swig arms during data recording. 

All signals were filtered and then sampled at 125Hz with 
12 bits precision per sample. Table I shows the specifications 
for all analog filters. 


TABLE I 

Specifications for analog filters. Li and fc 2 denote the lower 
AND UPPER cutoff FREQUENCIES AT 22DB ATTENUATION. 



ECG 

PPG 

Accelerometer 

Filter Type 

Band-pass 

Band-pass 

Low-pass 

/cl 

O.lHz 

0.25Hz 

- 

/c2 

lOOHz 

14.5Hz 

lOOHz 


In the simulation framework, ECG signals were only used 
as reference signals to extract the ground truth heart rate. One 
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channel PPG (site PPG(l) placed at the back of the wrist) 
and there-axis accelerometer were the data actually used. 
We downsampled (Nyquist sampling) PPG and accelerometer 
signals to 31.75Hz in the experiments. It should be noted that 
the accelerometer data contained aliasing when downsampled 
at 31.75Hz. However, such distortions were small as the rate 
of most activities during fitness training rarely exceed lOHz. 


B. Experiments Setup 

Fig. 8 shows the experiment setup. Quantized compressed 



remote-sensors receiver 


Fig. 8. The quantized compressed sensing was applied on PPG and three axial 
accelerometer data to reduce the transmission bit-budget. At the receiver, data 
were recovered from quantized measurements, and the TROIKA framework 
was applied on the recovered data to estimate the heart rate. 


sensing is applied on PPG and accelerometer signals. Raw 
segments {xi,X 2 ,X 3 ,X 4 } from each data channel are com¬ 
pressed by y = $x simultaneously. Sparse binary sensing 
matrices, whose entries consisted of only Os and Is, are used 
in the experiment. Such matrix preserves low-power property 
when implemented in FPGA [10]. We fix each column of $ 
consisting exactly 2 non-zero entries and also make sure $ is 
full row-rank in each iteration. The compressed measurements 
{yij y 2 , ya, Yt} are further quantized by Q(-) to reduce the 
transmission bit-budget. 

At the receiver, CS algorithms are used to recover signals 
from quantized measurements. The heart rate is estimated by 
the TROIKA [32] framework using the recovered PPG and 
accelerometer datasets. 

The codes and data reproducing the results in the experi¬ 
ments are available at https://github.com/liubenyuan/qsbl. 


C. Performance Measurement 

We use two performance metrics. One is the reconstruction 
SNR (RSNR), 


RSNR (dB) = lOlogio 


where x denotes the recovered signal of x. To assess the 
averaged performance over all segments, we refer to the 
average RSNR (ARSNR), defined by 

ARSNR(dB) = 101„g4|:(^^) ,37) 


where x^ is the ith signal segment and S is the total number 
of segments in a dataset. The second metric is the Structural 


SIMilarity index (SSIM) [33]. SSIM measures the similarity 
between recovered signals and original signals, which is a 
better metric than RSNR [15], [33]. 

The qualities of recovery are not only characterized by 
RSNR, but also by application specific requirements. There¬ 
fore, we perform a task-driven approach where the heart rate 
estimates from the recovered PPG and accelerometer signals 
are evaluated. The average absolute estimation error (Errorl), 
defined in [32], was, 

1 ^ 

Errorl = |BMPest(i) - BPMtrue(i)| (38) 

^ i=l 

where W is the total number of heart rate estimates, 
BPMest(*) is the estimated heart rate in the *th time window 
’ and BPMtrue(*) is the ground truth heart rate. The standard 
deviation of heart rate estimates, denoted by SDbpm, is 


SDbpm = 


\ 


Pearson correlation between the ground-truth and the heart rate 
estimates is also calculated for comparison. 


D. The Recovery Algorithms for Quantized CS 

Besides the proposed algorithm, we use the following two 
typical CS algorithms: QVMP [19] and BSBL-BO [3]. 

(1) QVMP [19]. QVMP is a variational Bayesian De- 
Quantization algorithm proposed in [19]. As shown in [19], 
it has better performance than QIHT [18] and LIRML [17]. 
However, QVMP can not recover physiological signals directly 
in time domain. Instead, we applied QVMP in transformed 
domain, where 

z = (€>’ip)0 -f e -f n 

is a sparse representation matrix and 6 is the sparse 
coefficients. QVMP firstly recovered 0, then x via x = \P0. 

In the experiment. Discrete Cosine Transform (DCT) matrix 
was selected for ’i' and we set A = 0.001/M, tol = le“®, 
maxiter = 400 for QVMP as it achieved best recovery 
performance on the datasets. 

(2) BSBL-BO [3], [14]. BSBL-BO [14] is the best perform¬ 
ing CS algorithm in recovering physiological signals such as 
fetal ECG [3] or EEG [15]. In the experiment, we observed 
that BSBL-BO can also recover signals from quantized mea¬ 
surements, by modeling quantization errors using a zero mean 
Normal distribution with a larger variance Ag. 

Throughout the experiment, BSBL-BO worked in non- 
sparse mode and directly recovered signals in time domain. We 
selected blkLen = 32, LearnLambda = 2, maxiter = 64 and 
LearnType = 1 for BSBL-BO as this setting achieved best 
recovery performance. The recoveries of BSBL-BO using real¬ 
valued measurement were also calculated, where the parameter 
LearnLambda = 0 was used. 

*The TROIKA algorithm [32] operates in a sliding window manner. A 
time window of T seconds is sliding on the signals with incremental step 
S seconds, and the heart rate estimates are based on the samples collected 
within this time window. We use default parameters (T = 8s, S = 2s) for 
TROIKA in our experiments. 
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Original signal 


Compressed Measurements 


Reconstructed by BSBL-BO (real-valued) 




Reconstructed by BSBL-BO (quantized) 






Fig. 9. Performances of different CS algorithms on quantized measurements, (a) a raw PPG segment x was collected in 4s (A^ = 128). (b) real-valued 
compressed measurements (M = 64, CR = 0.50) were generated via y = and then quantized by Q(-) with 2 bits per measurement. The transmission 
bit-budget was only MB = 128 bits, (c) on real-valued measurements, BSBL-BO recovered original signals with RSNR = 12.046 dB and SSIM = 0.957. 
(d) QVMP recovered signals in DCT domain with RSNR = 5.121 dB and SSIM = 0.784. (e)(f) BSBL-BO and BDQ recovered signals directly in time 
domain with RSNR = 5.174 dB, SSIM = 0.799 and RSNR = 7.422 dB, SSIM = 0.876 respectively. 


E. Results 

I) An illustrative example: To better understand the quan¬ 
tized compressed sensing and the quality of signal recovery 
using different CS algorithms, an illustrative example was 
given in Fig. 9. 

1) Fig. 9 (a). A raw PPG segment x of size N = 128 was 
shown. In our experiments, signals were divided into fix¬ 
sized segments and each segment was normalized (i.e., 
x/||x||) before compression. 

2) Fig. 9 (b). A segment x was compressed via CS, 
y = $x, where $ was a sparse binary matrix with 
M — QA rows and the compression ratio CR = 0.50. 
The compressed measurements were further quantized 
by Q(-) with B = 2 bits per measurement. The reference 
voltage for the quantizer was I^e/ = 0.70max(y). The 
quantized measurements z = Q(y) have only 4 finite 
values and contain saturations. The total transmission 
bit-budget was only MB = 128 bits. 

3) Fig. 9 (c) shows the recovery results of BSBL-BO using 
real-valued measurements y. 

4) Fig. 9 (d)(e)(f) shows the recovery results of QVMP, 
BSBL-BO and BDQ using the quantized measurements 
z. QVMP and BDQ used the prior information of 
quantization cell width A to assist recovery. 

Using the same settings in Fig. 9, we presented ARSNR and 
SSIM of different types of signals on the dataset ‘subject-02’ 
in Fig. 10. The ARSNR was calculated on all 74 segments 
in this dataset. The results in Fig. 9 and Fig. 10 shows that 
the proposed algorithm, BDQ, yielded better signal recoveries 
from quantized measurements than QVMP and BSBL-BO. 

We also presented results in recovering PPG signals using 
multiple quantization bit-depth B G {2,3,4,6,8} in Fig. 
11. The proposed algorithm, BDQ, was in average 1.640dB 
superior to BSBL-BO and 3.547dB to QVMP. Both BDQ and 
BSBL-BO were superior to QVMP. 

For larger bit-depth B such as B G {6,8}, the variance 
of the quantization error is small and can be approximated 


subject-02, PPG and Accelerometer Data subject-02, PPG and Accelerometer Data 




Fig. 10. The ARSNR and SSIM of signal recoveries w.r.t different types 
of signals. The parameters N = 128, M = 64, B = 2 were used. The 
three-axial accelerometer data were denoted by ax (x-axis), ay (y-axis) and 
az (z-axis) respectively. 


subject-02 



B 

Fig. 11. The ARSNR in recovering PPG signals in dataset ‘subject-02’ with 
multiple quantization bit-depth B. 


by a Normal distribution. However, in Fig. 11, we observed 
performance gap between BSBL-BO (real-valued) and BSBL- 
BO (quantized) when B G {6, 8}. This was largely caused 
by the saturation errors of the quantizer, whose distributions 
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are unbounded and can not be approximated by a Normal 
distribution. In contrast, BDQ yielded similar recovery per¬ 
formance to BSBL-BO (real-valued) when B € {6,8}. The 
reasons were two-fold, one was that the regularization (32)- 
(33) on the correlation matrix in BDQ can better exploit the 
highly correlated structure in PPG signals than the empirical 
method used in BSBL-BO, the other was the learning rules 
(35) for quantization errors can account for mild saturations 
introduced by the quantizer. 

2) The trade-off between the compression of CS and the 
quantization bit-depth: Quantized CS, when applied to data 
compression for low-energy telemonitoring, is basically a two- 
stage compressor. Firstly, we compress a raw segment of 
N samples via compressed sensing to achieve a preliminary 
compression ratio CR = {N—M)/N. Then the measurements 
are efficiently encoded by quantization with only B bits per 
measurement to further reduce the transmission bit-budget. 

We present in Fig. 12 the results with varying number of 
measurements M G {32,64,96} and quantization bit-depth 
B G {2,3,4,6,8}. We also present results of traditional 


The SSIM and the total transmission bit-budget (MB) 



Fig. 12. The SSIM w.r.t different transmission bit-budget MB in recovering 
the PPG channel of dataset ‘subject-02’. We fixed N = 128 and varied M S 
{32, 64, 96}. BDQ was used to recover signals from quantized measurements 
with different bit-depth B. 

sample-based compression using BSBL-BO, where the bits 
per compressed measurement is equal to the bit-depth of the 
signals in the dataset, i.e., B = Bi = 12. Fig. 13 shows the 
absolute heart rate estimation error (Errorl) and SDbpm with 
respect to different transmission bit-budget. 

From the results in Fig. 12 and Fig. 13, we found that both 
the SSIM and Errorl are affected by the configurations of 
M and B even if under the same transmission bit-budget. In 
fact, CS provides random mixing of signals to compressed 
measurements, each compressed measurement preserves in¬ 
formation on original signals while the quantization process 
drops information. Therefore, if we were going to quantize 
the compressed measurements with small bit-depth B, the 
compression ratio CR can not be high. In Eig. 13, the 
configuration M = 64 had the minimal Errorl for MB < 512 
bits. It is worth noting that smaller values of MR is always 
preferable since it reduces total transmission bit-budget. Eor 
N = 128, Ri = 12 and the optimal configuration (point A 


Heart rate estimation Error and the transmission bit-budget MB 



Transmission Bit-budget (bits) 


Fig. 13. The heart rate estimation error (En'orl) and SDbpm w.r.t different 
transmission bit-budget MB. For each configuration (M,B), the diameter 
of a solid circle is proportional to the SDbpm metric of this configuration. 
Point A denotes the optimal configuration where M = 64 and B = 2. 


in Eigure. 13) of quantized CS M = 64, R = 2, the bit 
compression ratio CR;,, defined in (14), is 




We achieved 91.67% bit compression ratio and transmitted 
128 bits instead of 128 samples for telemonitoring. 

3) Heart rate estimates from recovered datasets: We now 
presented results on the whole datasets using the optimal 
configurations of quantized CS with N — 128, M = 64, CR = 
0.50 and R = 2. Eig. 14 shows the absolute error (Errorl) 
per subject. Table II lists absolute error (Errorl), Standard 
Deviation (SDbpm) Pearson correlation calculated over all 
subjects. By jointly using the optimal configuration for quan¬ 
tized CS and the proposed algorithm, we achieved Errorl = 
2.596 (BPM), SDbpm = 3.625 (BPM) and 0.9899 Pearson 
correlation, which is closely to the result on non-compressed 
datasets. The results in Table. II also shows that for signal 
recovery from quantized measurements, BDQ was superior to 
both QVMP and BSBL-BO. 

The Scatter Plot between the ground-truth heart rates and 
the estimates using the recovered datasets is shown in Pig. 
15. The fitted line for BSBL-BO (CR = 0.50, R = 2) and 
BDQ (CR = 0.50, R = 2) was ?/ = 0.999a: - 0.866, (R^ = 
0.960) and y = 1.002a: — 1.333, (R^ = 0.980) respectively, 
where x indicates the ground-truth heart rate value and y is the 
estimates from recovered data, B? is a measure for goodness 
of linear fit. 


V. Discussions 
A. The quantizer Q(-) 

The performance of signal recovery from quantized mea¬ 
surements clearly depends on the choice of the quantizer. Eor 
the uniform quantizer used in this paper, the dilemma is the 
choice of reference voltages 14e/- 

In our experiments, signal x was normalized via x/||x|j and 
then compressed by y = 4&x, the reference voltage Vref for 
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Fig. 14. The absolute heart rate estimation eiTor (Errorl) on each subject using non-compressed datasets (denoted by ‘un-compressed’) and the recovered 
datasets by different CS algorithms. The solid and the dashed line denote mean±standai'd deviation of EiTorl averaged over all subjects using non-compressed 
datasets. Using the optimal configuration of the quantized CS (A^ = 128, M = 64, CR = 0.50, R = 2) and the proposed BDQ algorithm, we achieved 
stable and closely performance to that of ‘un-compressed’. 


TABLE II 

The AVERAGED ERRORI AND THE PEARSON CORRELATION OVER ALL 12 SUBJECTS. 



Uncompressed 

Real-valued CS 

BSBL-BO 

Quantized CS (CR = 

QVMP BSBL-BO 

0.50, B = 2) 

BDQ 

Enorl (BPM) 

2.464 

3.137 

3.355 

3.179 

2.596 

SDbpm (BPM) 

3.554 

5.022 

5.561 

5.030 

3.625 

Pearson Con'elation r 

0.9902 

0.9810 

0.9747 

0.9798 

0.9899 





(a) BSBL-BO (CR = 0.50, B = 2) (b) BDQ (CR = 0.50, B = 2) 

Fig. 15. Scatter Plot between the ground-truth heart rate and the estimates 
using the recovered data by BSBL-BO (CR = 0.50, B = 2) and BDQ 
(CR = 0.50, B = 2). The Pearson correlation was 0.9798 and 0.9899 for 
BSBL-BO and BDQ respectively. The ‘Fitted Line’ is the linear fit of the 
estimates to ground-truth heart rate. The value, which is an estimates for 
goodness of linear fit, was 0.960 for BSBL-BO and 0.980 for BDQ. 

the quantizer Q(-) was set to V-e/ = 0.70max(y). It is worth 
noting that Vref = 0.70max(y) is not the optimal reference 
voltage for this datasets and we do not search for such optimal 
values to avoid overfitting. However one should take care that 
for smaller values of Vref there may be more saturations, while 
for larger values of Vref there may be more underflows. 

At the first sight, signal normalization x/||x|| and max(y) 
are not practical for implementing in hardware and also for 
low energy applications. Instead, this problem is solvable by 
instead fixing the reference voltage for an ADC and using 
an Automatic Gain Control (AGC) circuit to tune the scale 
of signals in between [—Vref, Vref]- The gains of AGC must 
be transmitted alongside with the compressed bits for signal 


recovery. 

B. Signal recovery directly in the time domain 

In literature most CS algorithms recover signals in a trans¬ 
form domain where signals x can be sparsely represented. 
By suitably choosing the transformation matrix, one can 
improve the quality of recovery. 

However, in practice, physiological signals recorded by 
wearable devices are usually contaminated by various strong 
artifacts, such as artifacts due to body motion and hardware 
issues [3], [23]. As a result, many of these signals are less 
sparse in many known transform domains. Seeking/designing 
an optimal transformation matrix for a specific kind of physio¬ 
logical signals may be difficult. In this situation, recovering the 
signals from transform domains is not effective. In [3], BSBL- 
BO was used to recover raw fetal ECG recordings directly in 
the time domain by exploiting temporal correlations of these 
recordings, revealing that exploiting correlations of signals 
is an alternative method to exploiting sparsity. An obvious 
advantage of this method is that it avoids the seeking of an op¬ 
timal transformation matrix for recovery. Our proposed BDQ 
algorithm also suggested the effectiveness of this method. 

C. Quantized CS as data encoders for more types of signals 

The CS and the quantizer are analogue to the Discrete 
Wavelet Transform (DWT) and the scalar quantization unit 
used in the JPEG standard [34]. However, CS can be low- 
energy implemented in EPGA [10] and the additional 2- 
bit quantizer barely consumes any resource using fix-point 
arithmetic [10]. The data compressor proposed in this work 
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can be used as a low-energy data compressor/encoder for 
potentially more types of signals such as audios or images. 

VI. Conclusions 

In this paper, we present an approach for low-energy wire¬ 
less telemonitoring using quantized compressed sensing. The 
contributions of this paper are two-fold. First, we propose a 
two-stage data compressor, where signals are compressed by 
CS with a compression ratio CR = 0.50 and then quantized 
with B = 2 bits per measurement. Second, to pursue better 
signal recoveries from quantized measurements, we develop a 
Bayesian de-quantization algorithm that can exploit both the 
model of quantization errors and the correlated structure within 
signals. Experiment results showed that, by jointly using the 
proposed data compressor and the recovery algorithm, we 
achieved 2.596 absolute heart rate estimation errors and 0.9899 
Pearson correlation on whole datasets, which was closely 
to the performance on non-compressed datasets. The results 
imply that we can effectively transmit N bits instead of N 
samples, which may revolution the way we compress data for 
low-energy wireless telemonitoring. 
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